Relaxation and Zeno effect in qubit measurements 
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We consider a qubit interacting with its environment and continuously monitored by a detector 
represented by a point contact. Bloch-type equations describing the entire system of the qubit, the 
environment and the detector are derived. Using these equations we evaluate the detector current 
and its noise spectrum in terms of the decoherence and relaxation rates of the qubit. Simple ex- 
pressions are obtained that show how these quantities can be accurately measured. We demonstrate 
that due to interaction with the environment, the measurement can never localize a qubit even for 
infinite decoherence rate. 
PACS: 73.50.-h, 73.23.-b, 03.65.X. 



An account of decoherence and relaxation in quantum 
evolution of a two-level system (qubit), interacting with 
an environment and a measurement device, has become 
a problem of crucial importance in quantum computing. 
Numerous publications have appeared on this subject 
dealing with interactions either with a measurement de- 
vice (detector) [1-3] or with the environment (a thermal 
bath) [4,5]. Generally, the simultaneous influence of an 
environment and a detector on a qubit is very important 
for understanding qubit measurements because the envi- 
ronment and the detector act on the qubit in different 
ways. For instance, the environment at zero tempera- 
ture relaxes the qubit to its ground state. As a result 
the qubit finally appears in a pure state, even though it 
was initially in a statistical mixture. On the other hand, 
the measurement device puts the qubit in a statistical 
mixture, even if it was initially in a pure state. 

One of the most striking measurement effects in which 
the role of relaxation has not been investigated is the so- 
called Zcno paradox [6] . It consists of total freezing of a 
qubit in the limit of continuous measurement. Usually, it 
is associated with the projection postulate in the theory 
of quantum measurements. Indeed, it follows from the 
Schrodinger equation that the probability of a quantum 
transition from an initially occupied state of a qubit is 
P(At) = a(At) 2 , where a is a factor which depends on 
the system [6] . If we assume that At is the measurement 
time which determines the timescale on which the system 
is projected into the initial state, then after N successive 
measurements the probability of finding the qubit in its 
initial state, at time t = NAt, isP(t) = [l-a(At) 2 ](*/ Ai >. 
Thus P(t) — ► 1 for At — > Q, N — > oo and £=const. In- 
cluding the environment into the Schrodinger equation 
for the entire system one would expect from the above 
arguments that the relaxation processes could only affect 
the coefficient a, but cannot destroy the qubit localiza- 
tion in the limit of At — > 0. 

This conclusion, however, is not correct. We demon- 
strate in this Letter that any weak relaxation delocal- 
izes the qubit even in the limit of continuous measure- 



ment. It is shown by using new Bloch-type quantum rate 
equations for the description of a qubit interacting with 
a detector and its environment. These rate equations 
are derived from the microscopic Schrodinger equation 
for the entire system. Using these equations, we deter- 
mine the qubit behavior without any phenomenological 
parameters. In addition we are able to relate the qubit 
behavior to the detector outcome and therefore establish 
how decoherence and relaxation rates of the qubit can be 
accurately measured. 

With respect to the contradiction of our results to the 
Zeno paradox arguments, it is possible that the simple 
derivation described above is not valid when the envi- 
ronment is taken into account. One reason is that the 
expansion of P(t) = at 2 + ■ ■ ■ near t — is not valid be- 
cause of the discontinuity of P' it) at t = 0. The latter 
arises due to the transition to the continuous and suffi- 
ciently flat spectrum of the environment [7,8]. We present 
our result for an environment at zero temperature, where 
the interplay between decoherence and relaxation is most 
pronounced. (The detector, however is far from equilib- 
rium.) 

Let us consider the generic example of electrostatic 
qubit measurements. The qubit is an electron in a 
double-dot system, whereas the detector is a point con- 
tact placed near one of the dots, Fig. 1. When the elec- 
tron occupies the first well close to the point contact, 
Fig. la, the current is smaller than in Fig. lb due to 
the electron electrostatic repulsion. Thus, the electron 
is continuously monitored by the tunneling current. The 
entire system can be described by the tunneling Hamil- 
tonian H = Hq + Hpc + -Hint , where 

H = Eia\ai + E 2 ala 2 — floia\a 2 + a\a{) (1) 

is the qubit Hamiltonian and Hpc, Hi nt describe the 
point contact detector and its interaction with the qubit 

H PC = Eiajai + ^ E r a\a r + ^(Sl ;r (4a; + H.c.) 
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Here aj(ai) and a\.(a r ) are the creation (annihilation) 
operators in the left and the right reservoirs, and f2; r 
is the hopping amplitude between the states I and r of 
the reservoirs. For simplicity we consider electrons as 
spinlcss fcrmions and assume each of the reservoirs is at 
zero temperature. The interaction term Hi nt generates a 
change in the hopping amplitude, 6Qi r = Q'i r — Qi r . We 
assume that the hopping amplitude is weakly dependent 
on the states l,r, so that it can be replaced by its aver- 
age value, ~ n and 5£li r ~ <5f2. Thus the detector 
current is l\ = e2ir£l 2 PlPrV when the electron occupies 
the first dot, and I2 — e27r(f2 + Sfl) 2 PlPrV when the 
electron occupies the second dot [1]. Here pl,r are the 
density of states in the reservoirs and V = pl~ Pr is the 
bias voltage. 
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(a) (b) 

Fig. 1: A point-contact detector monitoring the elec- 
tron position in the double dot. pl,r denote the chemical 
potentials in the left and right reservoirs. 

It was shown in [1] that for the case of a large bias 
voltage V, one can reduce the Schrodinger equation for 
the entire system to Bloch-type rate equations describing 
the reduced density-matrix of the electron, aij(t). The 
diagonal terms of this density matrix, <J\i{t) and 022 (i) 
are the probabilities of finding the electron in the first 
dot or in the second dot, respectively. The off-diagonal 
matrix elements ("coherences"), <7i2(i), 021 (*) describe 
linear superpositions of these states. One finds [1] 

(Tn = -in (<7l2 - (721 ) (3a) 

(7i2 = ie<7i2 - iOo(2(7ii - 1) - (r d /2)cri2 , (3b) 



where 022 (i) = 1 — 011 (i) and 021 = f*2(i)- 
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Here 

E h2 = Te/2 and T d = U/h/e - \fhje} V/2ir is the 
decoherence rate due to interaction with the detector. 
Solving Eqs. (3a, 3b) one obtains that the decoherence 
term in Eq. (3b) leads to a vanishing of coherences in the 
limit t — + 00. Finally the electron density matrix becomes 
a completely random mixture for any initial condition: 

<Tij(t -> 00) = (l/2)5y. 

Now we introduce the environment, represented by a 
boson bath at zero temperature and interacting with the 
electron. First, consider the case in which the electron 
is not coupled to the detector, Hi nt — 0. Then, for 
any initial conditions, the electron relaxes to its low- 



est energy state (E + ) through boson emission. This 
final state can be found by disregarding any possible 
renormalization of the electron states due to interaction 
with the boson field. We thus diagonalize the Hamil- 
tonian Hq, Eq. (1) by using the following "rotation": 
a +! _ = ± cos(0/2)ai !2 +sin(#/2)a2,i, with tan# = 2£lo/e. 
Theni/o = {e/2)(ala--ala+), where e= (e 2 + 4f2g)i/2. 

In order to account for the above relaxation effects in 
the qubit evolution, we replace the Hamiltonian Hq by 
the following Hamiltonian: 

H' = H + Y, [E a bib a + c a {a\a-bl + aLa+b a )] , (4) 

a 

where b^ a (b a ) are creation (annihilation) operators of a 
boson with the energy E a . This Hamiltonian is essen- 
tially equivalent to the Lee model and has been investi- 
gated in numerous works [9]. In the weak coupling limit 
this model leads to the same results as the usual spin- 
boson model, although it includes an additional, direct 
coupling between the qubit states a\ 2 |0). 

If there is no interaction with the detector, Hi nt = 0, 
and c a is weakly dependent on E a , one can trace the bo- 
son variables in the Schrodinger equation, by reducing it 
to the rate equations for the electron density matrix in 
the basis states, aj_|0) [9] 

&—(t) = -r r cr__(t) (5a) 
(7+_ (i) = ^g(7 + _ (t) - (r r /2)<7+_ (t) , (5b) 

with <J ++ {t) = 1 - a (t) and <7_ + (t) = a\_(t). Here 

r r = 2 , Kc 2 Y p a is the relaxation rate and p a is the den- 
sity of boson states. Eqs. (5a, 5b) obviously reproduce 
an exponential decay (relaxation) of the electron from 
the upper level £L to the ground level E + : cr__(t) = 
cxp(— T r t). Note that the off-diagonal density-matrix el- 
ement <7-| (t) vanishes in the limit t — * 00, similar to 

(712 (t) in Eq. (3b). Yet, the disappearance of off-diagonal 
density matrix elements does not necessarily imply de- 
phasing. Indeed, in the case of relaxation, Eqs. (5a, 5b), 

the diagonal term a (t) vanishes as well for t — > 00. 

As a result, the qubit finally appears in a pure state (the 
ground state), in contrast to Eqs. (3a, 3b) leading to the 
statistical mixture. Now let us include the interaction 
with the detector. Then the Hamiltonian of the entire 
system becomes H' = H + Hpc + H int . It is useful to 
return to the original qubit basis a\ 2 |0), in which Hi nt , 
Eq. (2), has a simple form. The corresponding rate equa- 
tions for the qubit density matrix <7;j (t) are obtained by 
tracing the detector and boson degrees of freedom. These 
equations can be written directly by using the method of 
Refs. [1,10]. We obtain 

(Tn = -ifl (au - 021 ) - r r (Ke/2e)((7i2 + cr 2 i) 

- (T P /4)[1 + (e/e~) 2 ] (2(7n - 1) + r r (e/2e) (6a) 

(T12 = ie<7i2 - [iClo + r r (Ke/2e)] (2(7n - 1) 

+ T r [ K - (1/2)012 - K 2 ((7i2 + 021)] - {T d /2)a 12 , (6b) 
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where k = Qq/£. (Similar equations were obtained by 
Korotkov [11] in the weak coupling limit by using phe- 
nomenological arguments). 

Solving Eqs. (6a, 6b) for e = 0, one obtains the qubit 
density matrix for the stationary state (a = 0) 



a = a(t — > oo) 



1/2 
+ 2y) 



!//(! + 2») 
1/2 



(7) 



where y = T r /Td- This describes a heated qubit with 
an effective temperature T e ff — 2f2 /ln(l + Ay). This 
heating is caused by the measurement process [12]. 

Now, we investigate how the relaxation affects the 
time-dependence of the qubit density matrix. Consider 
again the symmetric case, e = 0. Let us evaluate the 
probability of finding the electron in the first dot, <7n(t). 
By solving Eqs. (6) for the initial conditions cth(O) = 1, 
a 12 (0) = 0. We find 



a ij( t ) = Y,n a ij (*)' wnere a ij = a ij ■ Starting from 
the microscopic Schrodinger equation for the entire sys- 
tem and using the same method as in Refs. [1,8,10] we 
can demonstrate that in the limit of high bias-voltage V 
of the detector, Fig. 1, the off-diagonal density-matrix el- 
ements oj™" ^ are decoupled from the diagonal elements 

er-"^ in the equation of motion [8]. As a result we arrive 
to the following Bloch-type rate equations: 
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in (a[f aft) + T r ( K e/2e)(a[f + t#>) 



+j03 1 -»S ) -w-7(«s , -''r') 
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-r r t/2 



• (™) • (™) ■<-> < (™) ( n )\ 



(C ie - 8 -' + C 2e - e + 4 ) (8) +r r [«C9'_^ + - (l/2)a™ - «Vi? + <)] 
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(9a) 



(9b) 



where e± = \(T d ± fi), ft = ^/T'j - 64ft 2 , and Ci, 2 = 
l±(r<j/ft). Thus, for the case of weak decoherence, Yd <C 
8fto, the electron displays damped oscillations between 
the dots with the Rabi frequency, \J (2ft ) 2 — (r^/4) 2 . 

For strong coupling to the detector, Yd S> 80o, the 
situation is different. If T r = 0, the electron would stay 
in the same dot for a long time ("quantum Zeno" effect). 
The dwell time t z obtained from Eq. (8) is T d /8ttQ. Thus 
the increase of T d leads to a freezing of the electron, which 
is totally localized in the limit of T d — > oo. This result 
is consistent with the Zeno paradox, based on the pro- 
jection postulate. Indeed, the "measurement time" At 
is inversely proportional to IV (See for instance [3,13]). 
Also one can observe from Eq. (8) that for small At, 
P(At) = [1 - o-ii (i)] oc (At) 2 . However, as follows from 
Eq. (8), the interaction with the environment essentially 
destroys the Zeno effect. We find that the Zeno time 
in this case is t z x = (IY/2) + (8fig/r d ). Therefore the 
continuous measurement cannot localize the electron for 
a long time, even if Td — > oo. The corresponding dwell 
time is restricted by the relaxation rate T" 1 . The disap- 
pearance of Zeno paradox can be understood as follows. 
It is crucial to observe that the relaxations due to the 
environment and the detector are factorizable, as seen in 
Eq. (8). As a result of this factorizability, an effective 
linear in At term, generated by purely exponential decay 
for small At, appears in the expansion of P(At). This 
leads to the elimination of the Zeno effect [7,8]. 

Now, we will relate the qubit behavior to the corre- 
sponding observable quantities. For this, we have to in- 
clude the detector states in the rate equations (6). We 

thus introduce the (reduced) density matrix uj"™ ' (t) , 
where index n denotes the number of electrons that 
have arrived in the right reservoir by the time, t [8]. 
This density matrix is related to the previous one by 



h+l. 



2e 



°12 + 



2 (n-1) 



(9c) 



where (3± = 1 ± (e/e) and 0'± = 1 ± (e/2e). Tracing 
Eqs. (9) over n and using o\\ it) + <J2(t) — 1 we obtain 
Eqs. (6). 

Eqs. (9) allow us to evaluate the average detector cur- 
rent and its shot-noise power spectrum. The (ensemble) 
average current is given by 



I(t) = e^nP n {t) = (h - h)(ru{t) + I 2 



(10) 



where P n (t) = (t) + (t) is the probability of find- 
ing n electrons in the collector by time t. As expected, 
the average current is directly related to the occupation 
of the first dot. The shot-noise power spectrum can be 
calculated via the McDonald formula [14,15] 

S(uj) = — [™ dtsm(wt)^N%(t) , (11) 
7r J Q at 

where N 2 R (t)=j: n n 2 P n (t). 

Now, we investigate how the relaxation and decoher- 
ence rates can be extracted from I(t) and S(io). Con- 
sider first the stationary detector current J = J(t — > oo) 
for the symmetric case (e — 0). It follows from Eqs. (7) 
and (10) that I = (I\ + h)/?, so that it is independent 
of y = Y r /Td- Therefore this ratio cannot be extracted 
from I for e = 0. However, for a non symmetric qubit, 
e 3> fto, the detector current becomes sensitive to y. In- 
deed, I=(h+ h)/2 for r r = and T d ^ , but I ~ h 
for Td = and T r ^ (since the relaxation puts the sys- 
tem into the lowest energy state, E + <~ E\ for e ^> fto). 
Using Eqs. (6) and (10) we obtain for e » ft 



I = AI 



y+(ft /e) 2 
y + 2{n /ef 



(12) 
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where AI = I\ — 1%. 

Although in the symmetric case e = 0, the relaxation 
does not affect the stationary current /. It instead affects 
its transient properties, which are reflected in the shot- 
noise spectrum of the detector current, S(u>), given by 
Eq. (11). One can write S(co) = So + AS(lu), where the 
first term So = e(I\ + 12) is the Schottky noise and the 
second term is the excess noise generated by the qubit 
dynamics. Generally the analytic expression for AS(uj) 
is rather lengthy. We therefore present it only for e = 
and in two limits: Td, T r <C ^0 and 1^ 3> f2o- ln the first 
case the excess noise can be very well approximated by 
a Lorentzian 



AS(u) 



(A/) 2 (r d + 2ry 



(T d + 2T r ) 2 + 16(w - 2f7 ) 2 



(13) 



This corresponds to the result of Korotkov, obtained by a 
Bayesian approach ( "continuous" wave function collapse) 
in the weak coupling limit [11]. 

The second case, Td 3> flo, corresponds to the Zeno 
effect regime. We find 



AS(lo) 



(Ai) 2 (mv d n 2 + r r r 2 T ) 



4[i> 2 + 4(^ 2 - 4tt 2 ) 2 ] + rvr d r 2 



2 ' 



(14) 



where T T = T d + T r 



is seen in Fig. 2, one can measure the relaxation rate of 
a qubit via the noise spectrum of the detector at zero 
frequency. Specifically, the relaxation rate can be lower 
by two or more orders of magnitude compared with the 
dephasing rate and still change the noise spectrum signif- 
icantly. With such sensitivity, one may more accurately 
measure the relaxation rate because this allows one to 
increase the output signal by increasing the coupling be- 
tween the qubit and the detector. 

In summary, we found that a qubit interacting with its 
environment and with a detector can be described by a 
set of modified Bloch-type equations in which the deco- 
herence and relaxation process are clearly distinguished. 
The most interesting result of our analysis is that there 
is no Zeno paradox when the relaxation due to the envi- 
ronment is taken into account. In addition, we obtained 
simple analytical expressions for the detector current and 
its noise spectrum. Using these findings, we proposed a 
new and possibly more accurate way to measure the qubit 
decoherence and relaxation rates. 
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ECS-0102500 and DMR-0121146. 



AS(Oi) 1.2 
(A/) 2 



WO) 




Fig. 2: The excess noise power spectrum of the detec- 
tor current for different decoherence and relaxation rates. 
The solid curves correspond to T r = and the dashed 
curves to T r = O.lfio- 

The excess noise, AS(w), for different values of r^IY 
and e — are shown in Fig. 2. As expected, for Td < Q , 
the Rabi oscillations generate a peak in the noise spec- 
trum at u> = 2f2o . The relaxation modifies this peak 
according to Eq. (13). In the case of large decoherence 
rate T d 3> ilo and T r = 0, the qubit is in the regime 
of the Zeno effect. This leads to a telegraph noise [11], 
resulting in a peak at w = 0. This peak, however, is 
strongly diminished in the presence of relaxation, even 
for a small T r , as given by Eq.(14). 

This strong dependence provides a new way to measure 
relaxation rate of a quantum system in experiments. As 
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